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We extend the one pion exchange model at quark level to include the short distance contributions 
coming from rj, a, p and U) exchange. This formalism is applied to discuss the possible molecular 
states of DD* / DD* , BB* / BB* , DD* , BB* , the pseudoscalar- vector systems with C = B = 1 and 
C = —B = 1 respectively. The "5 function" term contribution and the S-D mixing effects have been 
taken into account. We find the conclusions reached after including the heavier mesons exchange 
are qualitatively the same as those in the one pion exchange model. The previous suggestion that 
1 ++ BB* /BB* molecule should exist, is confirmed in the one boson exchange model, whereas DD* 
bound state should not exist. The DD* / DD* system can accomodate a 1 ++ molecule close to the 
threshold, the mixing between the molecule and the conventional charmonium has to be considered 
to identify this state with X(3872). For the BB* system, the pseudoscalar- vector systems with 
C = B = 1 and C = — B = 1, near threshold molecular states may exist. These bound states should 
be rather narrow, isospin is violated and the 1 = component is dominant. Experimental search 
channels for these states are suggested. 

PACS numbers: 12.39.Pn, 12.40. Yx, 13.75.Lb,12.39. Jh 



I. INTRODUCTION 



Since 1970s it is widely believed that Quantum Chromodynamics should accomodate a richer spectrum than just 
qq and qqq resonances, many possible nonconventional structures are suggested, e.g. glueballs (gg, ggg,...), hybrid 
mesons {qqg) and multiquark states (gggg, qqqqq, qqqqqq, qqqqqq). Unfortunately, so far there is still no uncontroversial 
evidence for nonconventional states experimentally except the hadronic molecules. The deuteron is a well-known 
example of hadronic molecule, and the approximate 10 5 known nuclear levels are all hadronic molecule. In the past 
few years, many new states have been reported, a striking feature is that some of them are close to the thresholds of 
certain two hadrons, which inspires the possible interpretation of hadronic molecule. 

Hadronic molecule is an old idea, about thirty years ago the possible hadronic molecules consisting of two charm 
mesons are suggested [lj], and -0(4040) was proposed to be a P wave D*D* molecule 0. Since in general molecule 
is weakly bound, the separation between the two hadrons in the molecule should be large. We can picture the two 
hadrons as interacting via a meson exchange potential At large distance, one pion exchange is dominant. Guided 
by the binding of deuteron, Tornqvist performed a systematic study of possible deuteronlike two-meson bound states 
[j, 0. The role of pion exchange in forming hadronic molecules was studied by Ericson and Karl Q . Recently Close et 
al. performed a pedagogic analysis of the overall sign, in addition they included the contribution of the "5 function" 
term which gives a 5 function in the effective potential when no rcgularization is used. In these original work, onl 
long distance one pion exchange has been considered, and the short distance contributions are neglected. In Rcf. 
Swanson assumed that the short distance dynamics is governed by the one gluon exchange induced constituent quark 
interchange mechanism, which results in state mixing. 

In the model of the nucleon-nucleon interaction, the long range part of the nucleon-nucleon force is quantitatively 
accounted for by the one tt exchange. However, the short and intermediate range interactions are governed by more 
complex dynamics. Combining the well-established one n exchange with the exchange of heavier bosons (e.g. scalar 
and vector mesons) to describe the behavior at short distance has been proved to be a very successful approach 
00, El- Physically, the scalar and vector meson exchange describes part of multiple pion exchange effects. For the 
two 7r exchange, if they interact and correlate in a P wave state, such a exchange can be modeled by p exchange. If 
the two correlated n pair is in a S-wave state, Durso et al. showed that one can approximate them by the exchange of 
a scalar a meson (l2| . Similarly, the correlated 3 7r exchange can be approximated by the exchange of one lu meson. 

Inspired by the nucleon-nucleon interaction, we shall represent the short distance interactions by the heavier bosons 
77, fj, p and to exchange instead of the quark interchange. The effective potential between two hadrons is obtained by 
summing over the interactions between light quarks or antiquarks as in the original work 0, [H, 0] . It is well-known 
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that one pion exchange between two light quarks results in two terms: the isospin dependent spin-spin interaction 
and tensor force. After taking into account the heavy bosons exchange, six additional terms appear including the 
spin-isospin independent central term, only isospin dependent term, isospin independent spin-spin interaction and 
tensor force, both isospin dependent and independent spin-orbit interactions. Consequently the situation becomes 
more complex than the only one pion exchange model. In our model, both the "<5 function" term and the S-D mixing 
effects would be considered, which have be shown to play an important role in the binding [H, @] ■ ^ n this work, we 
first give a good description of the deuteron in our model, which is an unambiguous hadronic molecule, then apply 
this formalism to the heavy flavor pseudoscalar-vector (PV) systems. Thus the predictions for the possible heavy 
flavor PV molecules are base on a solid and reliable foundation. This is a greater advantage over other approaches 
dealing with the dynamics of hadronic molecule, such as one boson exchange in the effective field theory [Tj, [l4[ and 
residual strong force with pairwise interactions [TEl . etc. 

The paper is organized as follows. In section II, the formalism of the one boson exchange model is presented, the 
effective potentials from pseudoscalar, scalar and vector meson exchange are given explicitly. In section III, we give 
the meson parameters involved in our model and the boson-quark couplings which are extracted from the boson- 
nucleon couplings. The formalism is applied to the deuteron in section IV, the DD* / DD* system and the molecular 
interpretation of X(3872) are investigated in section V. We further apply the one boson exchange approach to other 
heavy flavor PV systems in section VI, and possible molecular states arc discussed. Section VII is our conclusions 
and discussions section. The expressions for the matrix elements of the spin relevant operators are analytically given 
in the Appendix. 



II. THE FORMALISM OF ONE-BOSON EXCHANGE MODEL 



The construction of one-boson exchange interaction is constrained by the symmetry principle. To the leading order 
in the boson fields and their derivative, the effective interaction Lagrangian describing the coupling between the 
constituent quarks and the exchange boson fields is as follows 0, [l(| E| 

Pseudoscalar : C p = — g pqq ip(x)ij5ip(x)Lp(x) 

Scalar : C s = —g sqq if>(x)ip(x)(j)(x) 



Vector : C v = -g vqq ij(x)^(x)v»(x) - ^-j>(x)<T^(x)d»v v (x) 



Here m q is the constituent quark mass, ip{ x ) is the constituent quark Dirac spinor field, f{x), </>(x) and v^ix) are the 
isospin-singlet pseudoscalar, scalar and vector boson fields respectively. In this work we take m q = m u = ~ 313 
MeV, since we concentrate on constituent up and down quarks. If the isovector bosons are involved, the couplings 
enter in the form t ■ tp, t • cf> and r • respectively, where r is the well-known Pauli matrices. For the pseudoscalar, 

£ — 

another interaction term is allowed C' p = ^ 3 -ij){x)^ ,1 ^^'ij}{x)d ll Lp{x) , where m p is the exchange pseudoscalar mass, 
this Lagrangian has been used by Tornqvist [J, [f| and Close Q • By partial integration and using the equation of 
motion, one can easily show that C p and C' p are equivalent provided the coupling constants are related by 

fpqq _ 9pgg 
m p 2m q 

From the above effective interactions, the effective potential between two quarks in momentum space can be calcu- 
lated straightforwardly following the standard procedure. To the leading order in q 2 /m^, where q is the momentum 
transfer, the potentials are 

1 . Pseudoscalar boson exchange 

dlqq 0» ■ q)(<7j • q) 



4?7j2 q2 + ^2 

q 2 , q 2 ^(q) 



9pqq 



q 2 + q 2 + Mp 

,2 



(3) 



where Sy(q) = 3(0V - o)(&j • q) — en • &j, we have used fi p = m p — q$ instead of nip to approximately account 
for the recoil effect 0, H,0|- 
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2. Scalar boson exchange 



Vs(<l) 



algg iSi 3 ■ (P x q) 
2m 2 q 2 + fi 2 



(4) 



where Sy = |(<Xj + «7j), /i 
momentum. 

3. Vector boson exchange 
K(q) = 



9vqq 



9vqq 



q$ with m s the exchange scalar meson mass, and p denotes the total 



H~ ^Qvqqfvqq 



smj q 

Jvqqfvqq iSjj • (p X q) 



/' r 



2m 2 



q 2 + /i 2 



+ fvqq) 2 q 2 ^j(q) - 2q 2 (cr, • n 3 ) 



12m 2 q 
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(5) 



where fi 2 



<7o approximately reflects the recoil effect with m v the exchange vector meson mass. 



The effective potential in configuration space is obtained by Fourier transforming the momentum space potential. 

Kj(r)= (2^/ d3qelq '^ (q) 



(6) 



where i = p, s and v respectively. However, the resulting potentials are singular, which contains delta function, so the 
potentials have to be regularized. Considering the internal structure of the involved hadrons, one usually introduces 
form factor at each vertex. Here the form factor is taken as 



F(q) = 



A 2 -' 



A 2 - 



A ! 



x 2 



(7) 



where A is the so-called rcgularization parameter, m and q are the mass and the four momentum of the exchanged boson 
respectively with A" 2 = A 2 — q^. The form factor suppresses the contribution of high momentum, i.e. small distance. 
The presence of such a form factor is dictated by the extended structure of the hadrons. The parameter A, which 
governs the range of suppression, can be directly related to the hadron size that is approximately proportional to 1/A. 
However, since the question of hadron size is still very much open, the value of A is poorly known phcnomcnologically, 
and it is dependent on the model and application. In the nucleon-nucleon interaction, the A in the range 0.8-1.5 GeV 
has been used to fit the data. For the present application to heavy mesons, which have a much smaller size than 
nucleon, we would expect a larger regularization parameter A. We can straightforwardly obtain the effective potentials 
between two quarks in configuration space. For convenience, the following dimensionlcss functions are introduced. 



ex i f^i ' 

&xi 



Hq(A, m 

-ffl(A,?71 

H 2 {A 1 m ex ,fi 1 r 
H 3 (A,m ex ,fi,r 
Gi{A,m ex ,fi,r 
G 3 (A,m ex ,fi,r 



1 

fir 



1 

fir 



A 2 



L ex p—Xr 



1 



= 1 + — 



= 1+ — 



2 2 

fir fi~r z 



3 
fir 



1 



2fiX 
X(A 2 -mt x ) 
2fi 3 

1 ,X 1 



-Xr 



-Xr 



A 2 - mi 



-Xr 



fi 2 r 2 ' fir 



Xr fi fi 2 r 2 
3 

Xr + X 2 r 2 ' 



—e-v* — (1 H — - — I — ) 



1 

fi r 



cos(fir) - e Xr + 



X(A 2 



-V' 



r,3 



l exl c ~Xr 



r / ~ \ 3 sin (fir) 3cos{flr)-, 1 , 

- cos [fir) : — - (1 

ur fi z r z fir 



2fi 2 

X 2 e~ Xr 
fi 2 fir 



3 
X~r~ 



fir 

A 2 - 



2 f i 2 



X 2 



1 + Xr) 



-Xr 



X 2 r 2 fi 2 fir 



A 2 



2/""' 



r,2 



[l + Xr)' 



■,-Xr 



fir 



fir 



(8) 



Then the effective potentials between two quarks from one-boson exchange are 
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1 . Pscudoscalar boson exchange 



V P (r) 



2 -3 q 

Sri2^j[~ Gi(A,m p ,fi p ,r)tri ■ tTj + G 3 (A, m p , fx p , , (if = -p 2 p < 



(9) 



with S lj (r) = 3(<7j • f)(cr J • f) - er,; • CTj 
2. Scalar boson exchange 



K(r) 



Q 2 

47T 



H (A,m s ,n s ,r) 



8m 2 , 



Hi(A,m sl fi s ,r) 



2m\ 



H 2 (A,m sl p s ,r)L ■ Sj 



Here L = r x p is the angular momentum operator. 
3. Vector boson exchange 



V v (t) 



(JLv / 

47 V 



gLqHo(A,m v ,fj, v ,r) 



idvqq + ^9vqqfvqq)^v 



8m 2 q 



Hi(A,m v ,iJ, v ,r) 



{Svqq + /ti(J(j) 



H 3 (A,m v , fi v ,r)S l j{r) + 2Hi(A,m v , p Vl r)(cri ■ crj) 



-(3.9, 2 



gg 



qq > 12m q 

Ivqqfvqq) ^2 (A, Ttl v , p v , ?*)L • Sj.; 



(10) 



(11) 



For I ~ I isovector boson exchange, the above potential should be multiplied by the operator Ti-Tj in the isospin 
space. We have included the contribution of the "6 function" term in the above potentials, which gives the delta 
function when no regularization is used, since this contribution turns out to be important 0]. The effective 
potential between two hadrons are obtained by summing the interactions between light quarks or antiquarks 
via one boson exchange. 



III. MESON PARAMETERS AND COUPLING CONSTANTS 



As the well-known nuclear-nuclear interaction in the one boson exchange model, we shall take into account the 
contributions from pseudoscalar mesons tt and r\ exchange, that from scalar meson a exchange, and those from vector 
mesons p and lo exchange. The basic input parameters are the boson masses and the effective coupling constants 
between the exchanged bosons and the constituent quarks. The meson masses with their quantum numbers are taken 
from the compilation of the Particle Data Group [19J. For the constituent quark-meson coupling constants, one may 
derive suitable estimates from the phenomenologically known irNN, tjNN, <tNN, pNN and loNN coupling constants 
using the Goldbcrgcr-Treiman relation. Riska and Brown have demonstrated that the nucleon resonance transition 
couplings to tt. p and lo can be derived in the single-quark operator approximation (20| . which are in good agreement 
with the experimental data. Along the same way, we can straightforwardly derive the following relations between the 
boson-quark couplings and the boson-nucleon couplings, 



9nqq 


3 m q 

= 9-xnn, 

5 mjv 


gr/qq — 


m q 

g-nNN 

m N 






9pqq 


= gpNN, 


fpqq = 


3 m q 

7 JpNN 

5 tu n 




3 m q 

- 7 )9pNN 

5 m 7v 


9i*Jqq 


1 

— 2 g^NN, 
1 

— 2 g<jNN 


fuiqq = 


m q 

JljN N — 




m q 

— jg^NN 

m N 


9aqq 











where is the nucleon mass. In the present work, the constituent up(down) quark mass m u (d) is taken to be usual 
value — 313 MeV, which is about one third of the nucleon mass. The effective boson-nucleon coupling constants 

are taken from the well-known Bonn model [TTj ] . and a typical set of parameters is shown in Table [T] The uncertainty 
of the effective couplings will be taken into account later, all the coupling constants except g^NN would be reduced 
by a factor of two, since the experimental value for g^NN has been determined accurately from pion- nucleon and 



5 



nuclcon-nucleon scatterings. In the following, we shall explore the possible molecular states consisting of a pair heavy 
flavor pseudoscalar and vector mesons, their masses are taken from Particle Data Group [19| : m D o = 1864.84 MeV, 
m D ± = 1869.62 MeV, m D *o = 2006.97 MeV, m D *± = 2010.27 MeV, m B o = 5279.53 MeV, m B ± = 5279.15 MeV and 
m B , = 5325.1 MeV. 



Boson J G (J JJ ) 


Mass (MeV) 


3747T 




7T* l-(0") 


139.57 


14.9 




7T° l-(O-) 


134.98 


14.9 




j] + (0") 


547.85 


3.0 




a 0+(0+) 


600.0 


7.78 




p 1+(1") 


775.49 


0.95 


35.35 


w O-(l-) 


782.65 


20.0 


0.0 



TABLE I: Spin, parity, isospin, G-parity, the masses of the exchange bosons, and the meson-micleon coupling constants in the 
model. 



IV. DEUTERON FROM ONE BOSON EXCHANGE MODEL 



Deuteron is a uncontroversial proton-neutron bound state with J = 1 and 1 = 0. It has been established that 
the long distance one pion exchange is the main binding mechanism, and the tensor force plays a crucial role, which 
results in the 3 Si and 3 Di states mixing. Tornqvist and Close only considered the pion exchange contribution in 
Refs. 0, 0], however, the scalar meson a exchange and the vector mesons p, uj exchange turn out to be important in 
providing the short distance repulsion and the intermediate range attraction, consequently, we shall take into account 
the contributions from the heavier boson exchange in the following. The effective potential becomes 



V d (r) = V«(r) + V t ?(r) + V a d (r) + V«(r) + Vj(r) 

= VS(r) + V#(r)(<n • «r 2 ) + VftrXn • r 2 ) + V T rf (r)5 12 (f) + V&(r)(«7i • «r 2 )(Ti • r 2 ) 
+V$ I (r)S 12 (r)(r 1 ■ r 2 ) + V£ s (r)(L ■ S) + V L d s/ (r)(L • S)(n ■ r 2 ) 



(13) 



where S = \{cr\ + <r 2 ) is the total spin, and L is the relative angular momentum operator. In the isospin symmetry 
limit, the components V^(r), Vg(r) etc are given by 



Vg(r) 

V s d (r) 
Vf(r) 
VS(r) 

VSAr) 
V&(r) 

VtsAr) 



n 2 
47T 



Ho(A,m ai m ai r) 



8m 2 N 



Hi(A 1 m (T ,m a ,r) 



a 2 

47T 



in 
,3 



8m 2 N 



H 1 {A,m ul ,m ul ,r) 



9 n NN m r, v ( a \ 



(<), 



NN 



Unn) 2 ml 



in 



6m 2 N 



HA A 



9pNN rriK \ 9 2 P NN + 4 9pNNfpNN m 3 

. m p H Q {A, m p , m p , r) - — ZZZT H A A , m p, m P : r ) 



v -H 3 (A,m ri ,m r „r) - — 



in i2m 2 N 

2 3 
9tvNN m ir 

4-7T 12m^ 



HA^-,m w ,m„,r) 



8m 2 N 



NN + fuNN) 2 
(9pNN + fpNN) 2 TO p 

6m 2 



in 



in 



H 3 (A, m w ,m w ,r) 



■HA A ,m p ,m p ,', 



N 



tNN 



in i2m 2 N 



^ T H 3 (A,m 7r ,m^,r) 



(g P NN + fpNN) 2 m 

12m 2 H A^mp,m p ,r) 



in 



91nn <g 2 ( A; m ^ r) 3 &nn + WW ^g 2 ( A; m ^ r ) 



in 2m 2 N 
^9 2 p nn 



in 



pNN fpNN m 3 p 

2^? 



in 



■H 2 (A,m p ,m p ,r) 



2m 2 N 



(14) 



N 
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In the basis of 3 Si and 3 Di states, the deuteron potential can be written in the matrix form as 



V d = [vS(r) + V s d (r) 3Vf(r) W s d j(r)} ( J J ) + [^ SI (r) - 3V&(r)] 



'0 
1 



+ [vS(r)-W d I (r)} ^) (15) 

Taking into account the centrifugal barrier from D wave and solving the corresponding two channel Schrodinger 



equation numerically via the Fortran77 package FESSDE2.2 [211 ]. which can fastly and accurately solve the eigenvalue 
problem for systems of coupled Schrodinger equations, we find the binding energy Ed — 2.25 MeV for the cutoff 
parameter A = 808 MeV, and the corresponding wavefunction is presented in Fig. [TJ If we reduce half of the effective 
coupling constants except g n NN, the binding energy is found to be about 2.28 MeV with A = 970 MeV. From the 
wavefunction one can calculate the static properties of deuteron such as the root of mean square radius, the D wave 
probability, the magnetic moment and the quadrupole moment, which are in agreement with experimental data. We 
would like to note that the small binding energy of deuteron is a cancellation result of different contributions of opposite 
signs. The detailed results are listed in Table UH it is obvious the results are sensitive to the regularization parameter 
A, and the same conclusion has been drawn in the one pion exchange model 0, 0]- The binding energy variation with 
respect to A is shown in Fig. [2j the dependence is less sensitive than the one pion exchange model. It is obvious 
that the binding energy variation with A is dependent on the coupling constants. For the coupling constants listed in 
Table U the binding energy no longer monotonically increases with A in contrast with the one pion exchange model. 
To understand this peculiar behavior, we plot the three components of the deuteron effective potential in Eq. (|15|) 
in Fig. [3] We can see that both Vu(A, r) and V22(A, r) potentials arc repulsive, and they increase with A at short 
distance. However, at intermediate distance the relation |Vi2(A = 1.2GeV, r)| < |Vi2(A = 0.8GeV,r)| < |Vi2(A = 
0.9GcV, r)| < |Vi2(A = 1.6GcV, r)| is satisfied, the Vi2(A, r) doesn't monotonically increases with A. Therefore the 
non-monotonous behavior in Fig. [2ji mainly comes from the non-monotonous dependence of Vi2(A, r) potential on A, 
which is a cancellation result of various contributions. As has been discussed above, the heavy flavor system should 
admit a larger A than the deuteron. Therefore the above values of A with which the smaller deuteron binding energy 
is reproduced, would be assumed to be the lower bound in the following. 



A(MeV) 


e d (MeV) 


r rma (fm) P D : Ps(%) Md(MJv) 


Qd(fin a ) 


808 


2.25 


3.85 5.66:94.34 0.85 


0.27 


900 


5.33 


2.77 7.44:92.56 0.84 


0.20 


1000 


4.96 


2.87 7.37:92.63 0.84 


0.21 




all couplin; 


;s are reduced by half except g-^NN 




A(MeV) 


£d(MeV) 


rrms(fm) Pd : Ps Md(Miv) 


Qd(fm") 


970 


2.28 


3.84 6.52:93.48 0.84 


0.28 


1100 


5.65 


2.70 8.92:91.08 0.83 


0.20 


1200 


8.89 


2.28 10.26:89.74 0.82 


0.16 



TABLE II: The deuteron static properties in the one boson exchange potential model, where Ed is the binding energy, r rms 
is the root of mean square radius, Ps and Pd represent the S-state and D-state probabilities respectively, /id is the magnetic 
moment, and Qd denotes the quadrupole moment. 



V. POSSIBLE DD* /DD* HADRONIC MOLECULE AND X(3872) 

The narrow charmoniumlike state X(3872) was discovered by the Belle collaboration in the decay B + — > K + + 
X(3872) followed by X(3872) Jh/jn+ir- with a statistical significance of 10.3a [H]. The existence of X(3872) has 
been confirmed by CDF (23|, DO (2J| and Babar collaboration [25[. the CDF collaboration measured the X(3872) mass 
to be (3871.61 ± 0.16(stat) ± 0.19(sys.)) MeV. Its quantum number is strongly preferred to be 1 ++ (26|. In the one pion 
exchange model, Tornqvist suggested that X(3872) is a 1 ++ DD* / DD* molecule and isospin is strongly broken [la ]. 
Recently Close et al. re-analyzed X(3872) in the same model, the critical overall sign is corrected and the contribution 
of the "5 function" term is included Q. Swanson have taken into account both the long rang pion exchange and 
short range contribution arising from constituent quark interchange Q. Recently Zhu ct al. dynamically studied the 
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(a) (b) 

FIG. 2: The deuteron binding energy variation with respect to the regularization parameter A. Fig. [2K corresponds to the 
coupling constants shown in Table U and Fig. [2p for the couplings reduced by half. 



binding of X(3872) in the heavy quark effective theory [T3|. In this section, we will investigate the 1 ++ DD* / DD* 
system from the one boson exchange model at quark level, where the short range interactions are represented by the 
heavier bosons 77, a , p and to exchange instead of the quark interchange. 

There is only a sign difference (— 1) G between the quark-quark interaction and quark-antiquark interaction, and the 
magnitudes are the same, where G is the G-parity of the exchanged meson. The diagrams contributing to the DD* 
and DD* interactions are displayed in Fig. [4] Because of the parity conservation, DD* can only scatter into D*D 




012345 6 7 012345670123456 

rfGeV -1 ) r(GeV~') rtGeV -1 ) 



(a) (b) (c) 

FIG. 3: The three components of the deuteron effective potential in Eq. (I13|l , Fig. [3] a, Fig. [3] b and Fig. [jj] c respectively 
illustrate the Vn(A,r), Vi2(A, r) and V"22(A, r) components. The solid line, dashed, dotted and dash-dotted lines correspond to 
A = 0.8 GeV, 0.9 GeV, 1.2 GeV and 1.6 GeV respectively. 
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via the pseudoscalar n and ij exchange, and DD* scatters into DD* with the scalar a exchange, whereas the vector 
mesons p and ui exchange contribute to both processes. The effective potential for the 1 ++ DD* / DD* system is 

V x (r) = - V x (r)+V x (r) + V x (r)+V x (r)-Vj(v) 

= V x (r) + Vj?(r)(tn •<r j ) + V x (r)(r t ■ Tj ) + V x (r)S tJ (r) + V£<p k ,r)(tn ■ v^n ■ r s ) 

+V x (p k , r)Sy (f)fa ■ Tj ) + V x s (r)(L ■ Sy) + V x SI (r)(L ■ S«)fa ■ Tj ) (16) 



with i and j is the index of light quark or antiquark, p k (k = 1,2,3,4) takes four different values due to the mass 



difference within the D, D* and n isospin multiplets. The eight components V x (r), V x {r) etc are given by 



-,2 



An 

Vbjqq + ^9i*Jqqfu>qq TTl^ 

An 8m 2 



8m\' 



,2 



V x (r) = --^-m„ H (A, m a ,m a ,r) + -^^H 1 (A,m <J ,m (7 ,r) - ^^-m^H^A, m u , m u , r) 



An 



Hi(A,m u ,m u ,r) 



irX/ \ 9rjqq C5 „ / . \ (,9wqq fuiqq) ^7 tj / a \ 

47r Ylm z q 47r om~ 



II ( A ^ 9pqq + ^9pqqfpqq m p , 

^7 M = -^m p Ho(A,m p ,m p ,r) - — — J ^H x {A,m p ,m p ,r) 

An An 8mi 



drjqq rr / a \ . ^\9ujqq 4" fujqq) A f 7 

W = , 10 2 #3 A,m,,,^ 5 ,r) 4- ^- — 77- —H 3 (A, m u , /j. 7 , r) 

An 12m* An 12m i q 

-AiTi (A, m ff± ,o, M, r) - ^+W 2 g^ff^A, m p , M6 , r), p 2 > 



X I \ I 4jt 12mj 



V j(p r) 

1 g^L^G^A, m ff± ,o, A, r) - ^a^il! jg, g l(A , m p , Me , r), ,i 2 = - A 7 2 < 

^ (/ir) = 1 -#T^g 3 (A, m n ±,o,p,, r) - (^+/^) 2 ^ g 3 (A, m Pi M 6, r), M 2 > 

-4* lfej G 3(A, n^i.o.A.r) - (^+^) 2 i2 ^ g 3 (A, mp , M6 , r ), ^ = _~2 < 

u l / \ 5a 99 ™ct u / A \ . 3g^ qq + Ag^qqfuqq 

VlstM = — 4^2^ ''''^ 4- — ■^—^H 2 {A,m UJ ,m UJl r) 



^9pqq 4~ ^9pqqfpqq m p 

4-7T 2m 2 " 



^ /(r) = - ^ -JLg 2 (A, Wp , Wp , r) (17) 



where 



/x x = m 2 - (m_D«o - m D o) 2 , \i\ = m 2 ± - (m B »o - m D ±) 2 

111 \2 2 2 / \2 

2 2/ \2 2 2 / \2 

= m q - (m]j*o - m D o) , r J -6= m p- {m D *o - m D o ) 

Ms = m w - K" - ™r»o) 2 ( 18 ) 

These /Lt 2 parameters approximately represent the recoil effect due to different values of mu and mo* as in Rcfs. 
@j 0]- For the 77, er, p and w exchange processes, the mass difference of mo° and mu± as well as mn*° and mu»± 
are neglected, since they are much smaller comparing with m^, m p and m u . X(3872) is very close to the D°D*° 
threshold, however, it is about 8.3 MeV below the D + D*~ threshold. Hence, isospin symmetry is drastically broken 
fl7l. [l8j. For the J PC = 1 ++ DD* / DD* system, they can be in S wave or D wave similar to the deuteron, then the 
wavefunction of this system is written as 

|X(3872)) = ^1-L\(D°D*° + D°D*°) S ) + ^1-L\(D°D* + D°D*°) D ) 
r V2 r V2 

+ «3(r) 1 | (D+jD *- + D - D *+ )s) + ^M_L KD+D *- + (19) 
r v2 r V2 

where the subscript S and I? denote the system in S wave and D wave respectively. Ui(r), 1*2 (r), «3(r) and U4(r) are the 
spatial wavcfunctions. There arc four channels coupled with each other as has been shown above, and we might as well 
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FIG. 4: DD* and DD* interaction in one boson exchange model at quark level, where the thick line represents heavy quark or 
antiquark, and the thin line denotes light quark or antiquark. 



choose the basis to be II) 



A 1 \(D°D*° 



+ D°D*°) S ), |2) 



^.\(D°D*° 



D°D*°) D ), |3) = -L\{D+D*~ +D-D 



" 1 " " mh \"i — " 1 " ^ )uh\»i — " < " " )S1 

and 1 4) = -j=\{D + D*~ + D~ D* + )d) ■ Using the analytical formula for the matrix elements presented in the appendix, 
the effective potential for 1 ++ DD* / DD* can be written in the matrix form as 



V X {r) = [Vg(r) + V s x (r)] 



10 
10 
10 
1 



+V*(r) 



+V*> (r) I ° " 3 / 2 ° 

+v Lsin ooo 



to r)] 




LSI 



(r) 



-10-2 
0-10-2 
-2 0-10 
0-2 0-1 

/ y/2 2V2\ 

V2 -1 2^2 -2 

2\/2 V2 

V 2^2 -2 V2 -1 J 

/ 

3/2 3 



\ 3 3/2 



(20) 



In the above equation, the value of /zj? is fif for the up-left 2x2 matrix elements, and it is equal to /i| for the 
down-right 2x2 matrix elements. There is ambiguity in choosing /if, value for the processes D°D*° — ► D* + D~ or 
D + D*~ — > D*°D°, accordingly fif, can take the value /x| or n\ for the off-diagonal 2x2 matrix elements, the numerical 
results for both choices would be given in the following. The different fik values is due to the isospin symmetry breaking 
from the mass difference within the D, D* and 7r isospin multiplets. Taking into account the centrifugal barrier from 
D wave and solving the four channel coupled Schrodingcr equation using the package FESSDE2.2, the numerical 
results are listed in Table [TTTT It is remarkable that the 1 ++ DD* / DD* system could accomodate a molecular state 
with mass about 3871.6 MeV for A = 808 MeV, it is very close to the central value of_X(3872) mass 3871.61 MeV. 
The corresponding wavefunction is shown in Fig. O it is obvious that the D°D*° + D°D*° component dominates 
over the D + D*~ + D~ D* + component. Since the spatial wavefunctions u\(r) and u^{r) have the same sign, the same 
is true for 1*2 ( J") and 1*4 (r), thus the 7 = component in this state is predominant, it would be a isospin singlet in 
the isospin symmetry limit. From the results in Table IIII1 we notice that the predictions about the static properties 
for the two ji 2 choices are very similar to each other, and the difference is small. The isospin symmetry is strongly 
broken especially for the states near the threshold. The uncertainties induced by the effective coupling constants 
are considered, we reduce half of the couplings except g^NN, and the numerical results are given in Table HTT1 as 
well. For both choices of the coupling constants, the binding energy and other static properties are sensitive to the 
regularization parameter A, and the bound state mass dependence on A is displayed in Fig. [5J It is obvious that the 
bound state mass decreases monotonically with the regularization parameter A as in the one pion exchange model. In 
short summary, the predictions arc qualitatively the same as those in the one pion exchange model, even after we have 
included the contributions from 77, er, p and w exchange. Since unexpectedly large branch ratio of X(3872) — > t/j(2S)j 
recently was reported (27[, we have to take into account the mixing between the 1 ++ DD* / DD* molecule and the 
conventional charmonium state in order to identify this state with X(3872). This is outside the range of the present 
work. 
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M 2 


A(MeV) 


M(MeV) 


r rms (fm) 


Pif :PS J :P s t- :P£-(%) 




808 


3871.6 


7.02 


90.76:0.56:8.11:0.56 




840 


3870.4 


2.84 


78.23:1.08:19.59:1.11 


M§ 


850 


3869.8 


2.45 


75.26:1.21:22.29:1.23 




900 


3865.9 


1.61 


65.17:1.89:31.06:1.88 




1000 


3849.2 


1.08 


53.06:4.72:37.65:4.57 




808 


3871.7 


11.34 


94.40:0.38:4.86:0.36 




840 


3870.7 


3.19 


80 74-0 QQ-17 9fv1 m 


2 


850 


3870.2 


2.68 


77.44:1.14:20.28:1.15 




900 


3866.4 


1.66 


66.23:1.85:30.09:1.83 




1000 


3849.9 


1.09 


53.35:4.69:37.43:4.53 




all couplings except g n NN 


£ir6 r6duc6(l by lietlf 




A(MeV) 


M(MeV) 


r rms (fm) 






970 


3869.1 


2.13 


70.65:1.65:26.02:1.69 


,,2 
M3 


1100 


3860.1 


1.25 


0/1. OQ-Qfi QQ-O QZ. 

o / .z4.z.yo.oD.oo.z.yo 




1200 


3848.2 


1.00 


51.80:4.40:39.46:4.33 




970 


3869.5 


2.28 


72.56:1.57:24.28:1.60 


2 
M2 


1100 


3860.8 


1.27 


57.76:2.94:36.38:2.92 




1200 


3849.0 


1.01 


52.04:4.38:39.29:4.30 



TABLE III: The predictions about the mass, the root of mean square radius(rms) and the probabilities of the different compo- 
nents for the 1 ++ DD* / DD* molecule. 




FIG. 5: The four components spatial wavefunctions of the 1 ++ DD* / DD* system with A = 808 MeV. 



VI. POSSIBLE MOLECULAR STATES OF OTHER HEAVY FLAVOR PV SYSTEMS 

A. BB* jBB* system 

For the 1 ++ BB* jBB* system, the kinetic energy is greatly reduced due to the heavier mass of B meson, and the 
interaction potential has features similar to those of the DD* / DD* system except that the former is deeper than the 
latter. Therefore molecular states should be more easily formed. Following the same procedure as the DD* / DD* 
case, the numerical results are shown in Table HVl where the /i 2 ambiguity is considered. For the same value of A, 
the BB* /BB* system is really more strongly bound than the DD*/DD* system, its binding energy is a few tens of 
MeV, and the same was predicted in the one pion exchange model 0] and in the models [lj, [28| . It is obvious that 
the predictions about the static properties for the two /i 2 choices are approximately the same. We notice that the 
isospin symmetry breaking is less stronger than the charm system, this is because the mass difference of B° and B + 
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A (GeV) A (GeV) 



(a) (b) 

FIG. 6: The variation of the 1 ++ DD* / DD* bound state mass with respect to A. (a) corresponds to the coupling constants 
shown in Table [I] and (b) for the couplings reduced by half. 



is smaller than that of D° and D + as well as D*° and D* + . It is notable that there may be two molecular states for 
appropriate values of A. The corresponding wavefunctions for A = 1000 MeV and /j 2 = fil bl = m 2 ± — (mj. — mgo) 2 
are displayed in Fig. the first state is tightly bound, whereas the second is loosely bound. We notice that the first 
state is almost an isospin singlet, and the I = component is dominant for the second state. This state can no longer 
be produced through B meson decay because of its large mass, and we have to resort to hadron collider. We can 
search for this state at Tevatron via pp — > 7r + 7r~ Y(IS'), and LHC is more promising. 



B. DD* system with C = 2 



The interaction potentials arise from the one boson exchange between two antiquarks instead of a quark and 
antiquark pair, hence both the n exchange and u> exchange potentials have overall opposite sign relative to the 
DD* /DD* case. In this case wc have four coupled channels (D + D*°) s , (D+D*°) D , (D°D*+) S and (D°D*+) D . The 
numerical results are given in Table fVl For A = 808 MeV or A = 970 MeV, we find no bound state. A bound state 
with mass about 3873.9 MeV appears for A = 1600 MeV (about 3873.1 MeV for A = 1900 MeV if the couplings except 
9irNN are reduced half), and the corresponding wavefunction is shown in Fig. [5] We notice that the wavefunctions 
of 3 5i D + D*° and D a D* + have opposite signs, the same is true for the 3 Z?i D + D*° and D°D* + wavefunctions, 
therefore this state would be a isospin singlet in the isospin symmetry limit. We notice that the 3 Di probability are 
much larger than 3 Si probability for the state with A = 1800 MeV, although there is centrifugal barrier for the D 
wave state. Thus the S-D mixing effect induced by the tensor force is especially crucial for this state. In short, the 
bound state of the DD* system appears only for the regularization parameter A as large as 1600 MeV or 1900 MeV, 
which is beyond the range of 0.8 to 1.5 GeV favored by the nucleon-nucleon interaction. Moreover, the parameters 
that allow X(3872) to emerge as a DD* / DD* molecule exclude the DD* bound state, as can be seen from the results 
in section V. Consequently we tend to conclude that the DD* molecular state may not exist. 



C. BB* system with 5 = 2 



The situation is very similar to the DD* system except the different mass of D mesons and B mesons, we list the 
numerical results in Table fVTl Wc find a marginally bound state with mass 10603.9 MeV for A = 808 MeV, which is 
very close to the BB* threshold. Its binding energy is much smaller than that of the 1 ++ BB* /BB* system, however, 
the binding energy is less sensitive to A than the latter case. Fig. [5] displays the wavefunction of the bound state 
solution with mass 10602.3 MeV and A = 900 MeV. It is obvious that the B+B*° and B°B*+ wavefunctions have 
the opposite sign, then the 1 = component is dominant in this state. If the couplings except g^NN are reduced by 
half, a weakly bound state with mass about 10601.5 MeV is found as well assuming A = 970 MeV. These indicates a 
weakly bound BB* should exist, This is consistent with the results of Manohar and Wise form heavy quark effective 
theory [2{|. For a loosely bound molecule, the leading source of decay is dissociation, to a good approximation the 
dissociation will proceed via the free space decay of the constituent mesons. The spin-parity forbids its decay into 
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BB, therefore the BB* molecule is a very narrow state and it mainly decays into BBj. 

D. Pseudoscalar- vector system with C = B = 1 

This system could have the same quantum as B c meson or its antiparticle, and it is different from all the systems 
discussed above, eight channels instead of four channels arc coupled with each other under the one boson exchange 
interaction, i.e. {D+B*°) s , (D+B*°) D , (D°B*+) S , {D°B*+) D , (D*+B°) s , (D*+B°) D , (D*°B+) S and (D*°B+) D . 
We can investigate the possible bound states along the same line, although it is somewhat lengthy and tedious. 
There is ambiguity in choosing the p? value as well, for the DB* — > D*B scattering process, we could take p 2 = 
m ex ~~ ( m D' — mo) 2 or /i 2 = m 2 x — (ms* — tub) 2 , where m ex is the mass of the exchanged boson. Specifically for 
D + B*° — ► D* + B° via n exchange, we can choose p 2 = m 2 — (m D ,+ — m D +) 2 or p 2 = m 2 a — (m s »o — m B a) 2 . This 
ambiguity has been taken into account in our analysis. The numerical results are given in Table [VlTl For A = 808 
MeV, we find no bound state. With the choice p 2 = m 2 x — (to_d* — mo) 2 , a bound state with mass 7189.7 MeV is 
found for A = 850 MeV, However, this solution disappears if one chooses p? = m 2 ex — (raj. — tub) 2 - Only when A is 
around 880 MeV, the bound state solutions can be found for both p? choices. The difference of the static properties 
for the two p? choices is relatively larger than that of the above systems considered, this is because of the larger 
difference between m^. — m^, — 140 MeV and m^* — ms — 45 MeV. We notice that the D°B* + component has 
the largest probability in the states, since the threshold of D a B* + is lower than that of D + B*°, D* + B° and D*°B + . 
The wavefunction of the state with mass about 7185.9 MeV and A = 900 MeV is shown in Fig. [TU1 it is obvious 
all the eight components of the spatial wavefunction have the same sign, consequently this state would be isospin 
singlet in the isospin symmetry limit. Similar pattern of bound state solutions is predicted if the coupling constants 
except g^NN are reduced by half. This state is difficult to be produced, since both c and b have to be produced 
simultaneously. The direct production of this state at hadron collider such as LHC and Tevatron is most promising, 
and the indirect production via top quark decay is a possible alternative. Once produced, it should be very stable, 
DBtt and DBj are the main decay channels. 

E. Pseudoscalar-vector system with C = —B — 1 

The effective interaction potentials are induced by one boson exchange between two antiquarks, therefore both the 
7r exchange and u> exchange contributions give opposite sign between the C = B = 1 system and the C = —B = 1 
system, nevertheless the overall signs of r], a and p exchange potentials remain. We have eight coupled channels as well, 
(D+B*-)s, (D+B*-) D , {D°B*°) S (D°B *°) D , (D*+B-) s , (D* + B~) D , (D*°B°) S and (D*°B°) D are involved. The 
numerical results are given in Table IVlTTl It is remarkable that the p? and A dependence of the bound state solutions 
is similar to the C = B = 1 case. With the same p 2 and A values, the predictions for the static properties of the two 
systems are not drastically different from each other. Concretely for A = 900 MeV and p 2 = m 2 x — (mp* — m_o) 2 , 
we find a bound state with mass 7187.6 MeV for the C = —B = 1 system, and the mass of C = B = 1 bound state is 
7185.9 MeV, the difference is about 1.7 MeV. The corresponding wavefunction with A = 900 MeV is plotted in Fig. 
II H which can be roughly obtained by reversing the overall sign of the third, fourth, seventh and eighth components 
of the C = B = 1 system wavefunction in Fig. [TDJ To understand the similarity of the predictions for the C = B = 1 
and C = — B = 1 system, we turn to the one tt exchange model, the effective potential comprises a spin-spin potential 
proportional to (<7j • <jj){ri ■ Tj) and a tensor potential proportional to Sij(r)(Ti ■ Tj), where the isospin matrix Tj and 
the spin matrix er^ only act on the light quarks. In the basis of the eight channels listed above, these two operators 
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can be written as 8 x 8 matrices 



(<Ti -(Tj){Ti 
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(21) 



For the C = B = 1 pseudoscalar-vector system, the corresponding matrix representations are obtained by replacing 2 
and 2-\/2 with -2 and —2y/2 respectively in Eq. (j2Tj) . It is obvious both operators contribute to only the off-diagonal 
4x4 matrix elements. As a result, the eigenvalues of the corresponding Schrodinger equation for the C = B = 1 
and C = —B = 1 cases are exactly the same, if the small mass difference within the isospin multiplcts is neglected, 
and the eingen-wavefunction of one system can be obtained from another by reversing the overall sign of the third, 
fourth, seventh and eighth components. Therefore the heavy bosons 77, a, p and to exchange contributes to effective 
potential, and the pion exchange contribution is still dominant. In short summary, even after including shorter 
distance contributions from rj, a, p and lu exchange, the results obtained are qualitatively the same as those in the 
one 7r exchange model. The same conclusion has been reached for all the system consider above. 



VII. CONCLUSIONS AND DISCUSSIONS 



Motivated by the nuclcon-nucleon interaction, we have represented the short range interaction by heavier mesons 77, 
<7, p and lo exchange. The effective potentials between two hadrons are obtained by summing the interactions between 
light quarks or antiquarks via one boson exchange. The potential becomes more complicated than that in the one 
pion exchange model, and there are six additional terms which are proportional to 1, 73 • Tj, er.j • crj, SV, (r), L • Sij 
and (L • Sij)(ri ■ Tj) respectively. 

We first apply the one boson exchange formalism to the deutcron, then generalize to DD* / DD* , BB* / BB* , 
DD*, BB*, PV systems with C = B = 1 and C = —B = 1. S-D mixing effects has been taken into account, and 
the uncertainties from the regularization parameter A and effective coupling constants are considered. We find the 
conclusions reached are qualitatively the same as those in the one pion exchange model. This implies that the long 
range n exchange effects dominate the physics of a weakly bound hadronic molecule, and we can safely use one pion 
exchange model to qualitatively discuss the binding of molecule candidates. Since the predictions for the binding- 
energy and other static properties are sensitive to the regularization parameter A and the effective couplings, we are 
not able to predict the binding energies very precisely. If the potential is so strong that binding energy is large enough, 
we would be quite confident that such bound state must exist. However, the exact binding energy will depend on 
the details of the regularization and the effective couplings involved. Our results indicate that the 1 ++ BB* /BB* 
molecule should exist, whereas DD* bound state doesn't exist. For A = 808 MeV (970 MeV), the binding energy, 
D wave probability and other static properties of deuteron are produced, meanwhile near threshold 1 ++ DD* / DD* 
molecule is predicted. To identify this state with X(3872), the mixing between this DD* / DD* molecule and the 
conventional charmonium state should be further considered to be consistent with the recent experimental data on 
A(3872) -> ip{2S)-f H3|. For the BB* system, the PV systems with C = B = 1 and C = -B = 1, near threshold 
molecular states may exist. Similar to the 1 ++ DD* /DD* molecule, these states should be rather stable, isospin is 
drastically broken, and the / = component is dominant. Direct production of the above doubly heavy states at 
Tevatron and LHC is the most promising way. We can search for the 1 ++ BB* /BB* molecule via, pp — * 7r + 7r _ T(lS') 
at Tevatron. The BB* bound state mainly decays into BBj if it really exists. The dominant decay channels of the 
heavy flavor PV bound state with C = B = 1 arc DBir and DB-f, and the possible heavy flavor PV bound state with 
C = —B = 1 mainly decays into DBir and DBj. 
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In our model, the involved parameters include the effective quark-boson couplings, the masses of the exchanged 
bosons and the hadrons inside the molecule. Therefore this model is quite general, it can be widely used to dynamically 
study the possible molecular candidates. We will further apply the one boson exchange model to baryon-antibaryon 
system etc, and compare the predictions with the recent experimental observations [301 ] - 
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APPENDIX A: THE MATRIX ELEMENTS OF THE SPIN RELEVANT OPERATORS 

For initial state consisting of two mesons A and B, with relative angular momentum L, total spin S and total 
angular momentum J, its wavefunction is written as 

\{AB)LS, JMj) = J2 (LM L ;SM s \JMj}\LM L )\SMs) 
M L , M s 

( 1/2 1/2 S A ) 

= SaS b S 13 S 2 A 1/2 1/2 S B \\L(S 13 S 24 )S,JMj) (Al) 

Sl3,S , 24 I Sl3 S24 S 



where S = \f2S + 1. For the convenience of calculating the matrix elements of the spin-orbit operator L • S24. we can 
recouple the state as 

r t <? t i f V2 1/2 s A 

\(AB)LS,JMj)= (-V L+S+J SaSbS 13 S 24 SJ ls \ <T *f J Jf \l 1/2 1/2 S B 

x\(LS 2i )J LS S 13 ,JMj) (A2) 

In the same way, we can recouple the the final state \ (A'B')L'S', J'M'j) via the Wigner 6-j and 9-j coefficients. In the 
following, we shall present the matrix elements of four light quark operators involved in the work, which is helpful to 
calculating the matrix representation of the effective interactions. 

1. The unit operator 1 

Using Eq. (|Al[) . it is obvious that 



((A'B')L'S',J'M'j\1\(AB)LS,JMj) = 5 ll ,6 ss ,6jj,5 Mj m> J2 S A S' A S B S' B S 2 13 Sl 4 

£'i3,S24 

r 1/2 1/2 s A ) { 1/2 1/2 s' A ) 

x I 1/2 1/2 S B \ ! 1/2 1/2 S' B \ = 5 ll ,8 ss ,5 SaS>a 5 SbS , b 5jj,5 MjM , (A3) 

I 5*13 5*24 S I 5*13 5*24 S 



2. The spin-spin operator er 2 ■ cr 4 

((A'B')L'S',J'M'j\<t 2 -<t 4 \(AB)LS,JMj) = 8 ll ,6ss'Sjj>S Mj m> j S A S A SBS B Sf 3 S% 4 

1/2 1/2 S A ) ( 1/2 1/2 S' A 
x[25 24 (5 24 + l)-3] { 1/2 1/2 S B W 1/2 1/2 S' B }> (A4) 

5*13 $24 S J [ S*i3 S*24 S 

where the spin operators er 2 and er 4 only act on the light quarks and antiquarks 
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3. The spin-orbit operator L • S24 

((A'B>)L'S',J'M' j \L.S 2 4\{AB)LS,JMj} = 5 ll ,5jj,5 Mj m' j J2 (-l) s+s ' +2L+2J S A S' A 

xS b S' b SS'S 2 13 SIJ 2 ls ^[J ls (J ls + 1) - L(L + 1) - S 2i (S 2i + 1)] I S ' 24 Jl s s 
L So, Jr,l f 1/2 1/2 Sa ) \ 1/2 1/2 ^ 1 

5i 3 J S> 1/2 1/2 5s 1/2 1/2 5 * (A5) 

^ [ 5l3 S24 S J ^ 'S'24 5 J 

where S24 = \{&2 + ""4); L is the relative spatial angular momentum. The matrix elements of L • S24 can be 
calculated by the Wigner-Echart theorem 3l|, and the same result has been obtained. 

4. The tensor operator 524(f) = 3(cr 2 • r)(<X4 • f ) — <r 2 • er 4 

It can be checked that the tensor operator 524 (r) is proportional to the scalar product of two rank-2 tensor 

(2) 

operators Y 2m and S m with m = 0,±1,±2, where Y 2m is the spherical harmonic function of degree 2, and the 

( 2 ) 

five components of S m are 

S2 ■* = 2*52+54+, S^ ~ — ~ O^ 20 '^- + S2+S40), Sq 2 ^ = — ^-(52 - 54+ — 45 2 o54o + 52+ 54_) 

S^l — —(S2-S40 + S20S4-), S^ 2 \ = —S2-S4- (A-6) 

Here 52+ = \{p2x + i^y), S20 = \ a iQ and 52- = \{&2x — i&2y)- The same convention applies to S^± and 540, 
the spin operators er 2 and cr 4 only act on the light quark and antiquarks. Using the Wigner-Echart theorem, 
the matrix element of this tensor operator can be obtained, although it is somewhat lengthy. 

((A'B')L'S', J'M'j\S M (t)\(AB)LS, JMj) =5jj,5 MjM ^Vw ]T 6 Sa4 , 1 {-l) J+L+L '+* 8 '+ 3 » +8 » 

1/2 1/2 S A 

1/2 1/2 S B 

5i3 5 2 4 5 

1/2 1/2 S' A 

x <; 1/2 1/2 s' B } 

5i3 524 5' 

The above expression is apparently different from the results in Ref. Q, However, the numerical results of all 
the matrix elements arc the same. 
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(a) (b) 

FIG. 7: The spatial wavefunctions of 1 ++ BB*/BB* molecule with A = 1000 MeV. There are two bound states with mass 
10457.6 MeV and 10600.4 MeV respectively, Fig. [TJi is for the first state, and the Fig. [7p for the second. 




FIG. 8: The wavefunction for the DD* system assuming A = 1600 MeV and fi 2 = fi 2 cl = m^ — (m D ,+ — m D +) 2 , its mass 
approximately is 3873.9 MeV. 
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FIG. 9: The spatial wavefunction for the BB* system assuming A = 900 MeV and fj, 2 = fj, 2 bl = m^o — (wis* — m-B+) 2 , its mass 
is 10602.3 MeV. 



0.7 



v 0.4 
E 

"3" 0.3 





*S, D*B'° 




! D, D*B'° _ 




s s t d°b"* 




'd, d"b'* " 




"8 D'*B° 




"D, D'*B° 




'S, D"°B* - 




"D D'"B* 




. :: 





1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
r(fm) 



FIG. 10: The wavefunction of the C — B — 1 pseudoscalar-vector system with A = 900 MeV and fi 2 = m 2 x — (m_o* — mo) 2 , 
the mass of this state is about 7185.9 MeV. 




FIG. 11: The wavefunction of the C — — B — 1 pseudoscalar-vector system with A = 900 MeV and fi 2 = m 2 x — (run* — mn) 2 , 
its mass is about 7187.6 MeV. 
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r 


'+ BB* /BB* 


M 


A(MeV ) 


JVL(MeV ) 


r rms (fm) 


P^:PS J :P+-:P+-(%) 




808 


10565.3 


0.60 


47.70:2.05:48.20:2.05 


mL 


900 


10543.5 


0.59 


44.11:5.69:44.50:5.70 


1000 


10457.6 


0.52 


27.52:22.41:27.64:22.44 






10600.4 


1.73 


37.10:9.43:44.26:9.22 




808 


10565.3 


0.60 


47.70:2.05:48.20:2.05 


2 

^bb2 


900 


10543.5 


0.59 


44.11:5.69:44.49:5.70 


1000 


10457.6 


0.52 


27.52:22.41:27.64:22.44 






10600.4 


1.72 


37.10:9.43:44.25:9.22 




all coupling except g^NN are 


reduced by half 




A(MeV) 


M(MeV) 


r rms (fm) 


Vf : Pff : P+- :P+-(%) 




970 


10544.8 


0.55 


45.24:4.60:45.55:4.61 


*4l 


1100 


10503.9 


0.51 


40.12:9.78:40.30:9.80 


1200 


10443.1 


0.46 


33.66:16.27:33.77:16.29 






10601.9 


1.91 


38.82:8.03:45.26:7.89 




970 


10544.8 


0.55 


45.24:4.60:45.55:4.61 


2 

L l bb2 


1100 


10503.9 


0.51 


40.12:9.78:40.30:9.80 


1200 


10443.1 


0.46 


33.66:16.27:33.77:16.29 






10601.9 


1.91 


38.83:8.04:45.25:7.89 



TABLE IV: The predictions about the mass, the root of mean square radius(rms) and the probabilities of the different compo- 
nents for the 1 ++ BB* /BB* system with = m ^± ~ ( m s* — m B ) 2 an d Mbt2 ~ m ?r± — ( m s* — m s+) 2 - 



DD* system with C = 2 




A(MeV) 


M(MeV) 


r rms (fm) 


P+":P+ U :P^:P U D +(%) 




1600 


3873.9 


3.15 


34.61:4.03:56.82:4.54 


2 

Mccl 


1700 


3865.1 


1.30 


20.29:28.37:22.63:28.71 




1800 


3770.9 


0.58 


0.19:49.78:0.19:49.84 






3872.8 


2.48 


41.47:1.38:55.54:1.63 




1600 


3873.9 


3.16 


34.54:4.04:56.87:4.55 


2 

Mcc2 


1700 


3865.1 


1.30 


20.31:28.35:22.65:28.70 




1800 


3771.0 


0.58 


0.19:49.78:0.19:49.84 






3872.8 


2.49 


41.41:1.38:55.58:1.63 




all couplings are 


reduced by half except g^NN 


•2 

M 


A(MeV) 


M(MeV) 


r rms (fm) 


P+ U :Pr:P^:P U o + (%) 




1900 


3873.1 


2.53 


36.05:5.90:51.55:6.50 


2 

Mccl 


2000 


3870.1 


1.82 


38.09:8.08:45.32:8.52 




2100 


3865.3 


1.44 


37.67:10.21:41.57:10.55 




2200 


3858.4 


1.19 


36.31:12.42:38.60:12.68 




1900 


3873.1 


2.54 


36.01:5.91:51.58:6.51 


2 

Mcc2 


2000 


3870.1 


1.82 


38.07:8.08:45.32:8.53 




2100 


3865.4 


1.44 


37.66:10.22:41.57:10.55 




2200 


3858.4 


1.19 


36.30:12.42:38.59:12.68 



TABLE V: The predictions about the mass, rms and the probabilities of the different components for the DD* system, 
and P^"° denotes the probability of S wave D + D*° in the state, and the meaning of P+°, P° s + and P£+ is similar. Here 
Mcci = "Vo - (m D ,+ - m D +) 2 and /4 C 2 = "i^o - (m D ,a - m D o) 2 
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BB 


system with B = 2 




A(MeV) 


M(MeV) 


rims (&a) 


P+ u :P+ u :Pii+ : P^(%) 




808 


10603.9 


4.09 


59.37:6.74:28.49:5.40 




900 


10602.3 


2.23 


43.52:10.84:35.63:10.02 




1000 


10598.8 


1.61 


37.22:14.37:34.48:13.93 




1100 


10592.2 


1.27 


31.33:19.34:30.24:19.09 




808 


10603.9 


4.09 


59.38:6.74:28.48:5.40 


2 

t^bb'Z 


900 


10602.3 


2.23 


43.52:10.84:35.63:10.01 




1000 


10598.8 


1.61 


37.22:14.37:34.48:13.93 




1100 


10592.2 


1.27 


31.33:19.34:30.24:19.09 




all couplings except g n NN are 


reduced by half 




A(MeV) 


M(MeV) 


r rms (fin) 


P+ U :P+ U :P^:P^(%) 




970 


10601.5 


1.99 


39.93:13.08:34.67:12.31 


2 

^bbl 


1000 


10600.7 


1.82 


38.34:14.02:34.30:13.34 




1100 


10596.6 


1.44 


34.37:16.80:32.47:16.36 




1200 


10590.3 


1.19 


31.33:19.33:30.32:19.03 




970 


10601.5 


1.99 


39.94:13.08:34.67:12.31 


2 

I- L bb2 


1000 


10600.7 


1.82 


38.34:14.02:34.30:13.34 




1100 


10596.6 


1.44 


34.37:16.80:32.47:16.36 




1200 


10590.3 


1.19 


31.33:19.33:30.32:19.03 



TABLE VI: The predictions about the mass, rms and the probabilities of the different components for the BB* system with 
B=2. Pg° represents the probability of S-wave B + B*°. Here nl bl = mjj. — (tub* — m B +) 2 and /j% b2 = m^ — (ms* — m B o) 2 . 



The pscudoscalar-vector system with C — B — 1 




A(MeV) M(MoV) r rm3 (fm) P+°(DB*) : P+°(DB*) : P°+(DB*) : P° + (DB*) : P+°(D*B) : P + °(D*B) : P°+(D*B) : P°+(D*B)(%) 




850 7189.7 5.54 8.10:0.02:89.69:0.02:0.87:0.36:0.67:0.26 
880 7187.9 2.05 21.88:0.07:72.40:0.08:2.09:0.96:1.73:0.80 
900 7185.9 1.58 27.34:0.11:65.19:0.12:2.54:1.35:2.16:1.19 
1000 7157.2 0.86 36.55:0.01:46.27:0.01:1.39:7.32:1.21:7.26 


2 


850 no bounded 

880 7189.5 4.11 10.70:0.02:86.98:0.02:0.79:0.50:0.61:0.39 
900 7188.1 2.18 20.59:0.04:75.09:0.04:1.34:0.98:1.10:0.83 
1000 7161.1 0.88 36.86:0.22:47.18:0.23:0.48:7.33:0.41:7.31 


all couplings are reduced by half except g-^NN 




A(MoV) M(McV) r rm3 (fm) P+°(DB*) : P+°(DB*) : Pg+(DB*) : P°+(DB*) : P+°(D*B) : P + °(D*B) : P[!+(D*B) : P°+(D*B)(%) 


2 


970 7189.2 3.13 16.76:0.07 
1000 7187.6 1.88 25.51:0.13 
1100 7177.1 1.03 34.50:0.42 
1200 7156.3 0.78 34.62:0.72 


78.38:0.08:1.75:0.85:1.45:0.67 
66.52:0.15:2.78:1.35:2.42:1.15 
49.32:0.45:4.94:2.96:4.62:2.81 
41.51:0.75:5.85:5.51:5.64:5.42 


2 

^Sc2 


970 no bounded — 

1000 7189.6 4.90 11.19:0.03:86.02:0.03:0.90:0.62:0.74:0.49 
1100 7181.9 1.20 34.04:0.20:54.01:0.21:3.26:2.67:3.08:2.53 
1200 7163.0 0.83 36.21:0.33:43.97:0.34:4.08:5.57:4.00:5.51 



TABLE VII: The predictions for the static properties of the PV system with C = B = 1, where Pg °(DB*) denotes the probability 
of S-wave D + B*°, and Pg°(D*B) denotes the probability of S-wave D* + B°. Here /i 2 cl = m ex — (mr>* — mo) 2 and [i 2 c2 = 
m 1x — (rn b* — tub) 2 with m ex the exchanged boson mass. 
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The pscudoscalar- vector system with C — —B — 1 




A(MeV) M(McV) r rms (fm) P + ~(DB*) : P + ~(DB*) : P^DB*) : P™(DB*) : P + ~(D*B) : P + ~(D*B) : Pg°(D*B) : P°°(D*B)(%) 




880 7189.2 3.11 15.16:0.04:80.15:0.05:2.21:0.45:1.63:0.32 
900 7187.6 1.85 23.82:0.08:68.07:0.09:3.83:0.61:3.02:0.48 
1000 7169.0 0.79 30.56:0.28:46.30:0.30:11.86:0.74:9.27:0.69 
1050 7148.4 0.51 15.57:0.02:53.99:0.04:22.84:0.10:7.41:0.03 
7154.3 0.63 53.57:0.36:16.62:0.34:5.65:0.59:22.22:0.66 


2 


880 no bounded 

900 7189.7 5.17 9.30:0.02:88.09:0.02:1.25:0.27:0.87:0.19 
1000 7176.2 0.91 30.77:0.20:50.54:0.21:9.62:0.70:7.32:0.65 
1050 7154.6 0.52 25.22:0.00:45.64:0.01:17.96:0.04:11.14:0.00 
7163.3 0.70 46.16:0.31:27.72:0.30:7.95:0.64:16.25:0.68 


all couplings except g-^NN are reduced half 




A(McV) M(McV) r rms (fm) P + ~(DB*) : P + ~(DB*) : P f j! (DB*) : Pg>(DB*) : P + ~(D*B) : P + ~(D*B) : Pg°(D*B) : Pg 1 (D*B)(%) 


M 2 d 


970 7189.4 3.57 15.68:0.06:78.98:0.07:2.29:0.58:1.90:0.44 
1020 7185.6 1.38 30.23:0.20:56.36:0.21:5.82:1.03:5.25:0.90 
1100 7173.2 0.82 33.18:0.42:43.08:0.45:10.53:1.24:9.93:1.17 
1200 7147.4 0.59 30.96:0.68:35.69:0.70:15.05:1.28:14.40:1.25 


2 

Mi,e2 


970 no bounded 

1020 7189.2 3.03 18.64:0.07:75.25:0.07:2.60:0.63:2.24:0.51 
1100 7180.5 1.00 33.61:0.29:47.78:0.30:8.09:1.16:7.68:1.09 
1200 7158.2 0.64 32.34:0.58:37.41:0.60:13.43:1.28:13.09:1.26 



TABLE VIII: The predictions about the static properties of the PV system with C — —B — 1, where P,t~(DB*) denotes the 
probability of S wave D + B*~ , and Pg~(D*B) denotes the probability of S wave D* + B~ . Here (jl 2 c1 — m 2 x — (mo* — thd) 2 
and fj. 2 c2 = m ex ~ {ttib' — ttib) 2 with m ex the exchanged boson mass. 



